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The paper is concerned with the problem on rolling of a homogeneous ball on an arbitrary surface. New cases when the 
problem is solved by quadratures are presented. The paper also indicates a special case when an additional integral and 
invariant measure exist. Using this case, we obtain a nonholonomic generalization of the Jacobi problem for the inertial 
motion of a point on an ellipsoid. For a ball rolling, it is also shown that on an arbitrary cylinder in the gravity field 
the ball's motion is bounded and, on the average, it does not move downwards. All the results of the paper considerably 
expand the results obtained by E. Routh in XIX century. 
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1. Introduction 

In this paper we consider the problem of sliding-free rolling of a dynamically symmetrical ball (the 
central tensor of inertia is spherical I = /^E) on an arbitrary surface. As it was indicated by E. Routh 
in his famous treatise jlUj , if the surface is a surface of revolution the problem is integrable even in the 
presence of axisymmetric potential fields. Here we give a more complete analysis of the Routh solution 
for the body of revolution, and present new integrals for the case of ball rolling on non-symmetrical 
surfaces of the second order. 



2. The equations for a ball moving on a surface 



In rigid body dynamics it is customary to introduce a body-fixed frame of reference. However, while 
studying motion of a homogeneous ball, it is more convenient to write the equations of motion with 
respect to a certain fixed in space frame of reference. In such a frame a balance of linear momentum 
and of angular momentum with respect to the ball's center of mass, involving the reaction and the 
external forces, may be written as 



mv = N + F, 



a x N + M F . 



(2.1) 



The sliding-free condition reads (the contact point velocity vanishes) 




v + u x a = 0. 



(2.2) 



Fig. 1. The rolling of a ball on 
a surface (G is the center of 
mass, Q is a point of contact 
with the surface) 



Here m is the mass of the ball, v isthe velocity of its center of mass, u is 
the angular velocity, I = fiE is the (spherical) central tensor of inertia, 
a is the vector from the center of mass to the point of contact, R is 
the ball's radius, N is the reaction at the contact point (see Fig. F 
and Mp are the external force and the moment of forces with respect to 
the point of contact respectively. 

Eliminating from these equations the reaction N and making use of 
the fact that the contact point velocities on the surface and on the ball 
coincide, we obtain the system of six equations: 



M = Dj x (w x 7) + Mp, r + Rj = u x Rj. 



(2.3) 



Here D = mR 2 . These equations govern the behaviour of the vector of the angular momentum (with 
respect to the contact point M), and of the vector 7 = —R~ 1 a normal to the surface (Fig. Q). The 
vectors u and r (the position vector of the contact point) are obtained from the relations 



M = fxuj + D~{ x (w x 7), 7 



VF(r) 
\VF(r) 



(2.4) 
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F(r) = is the equation of the surface on which the ball rolls. The last equation in (|2.4|) defines the 
Gaussian map. Further on, following Routh, we will use the explicit form of the surface on which the 
ball's center of mass is moving. This surface, defined by the position vector r' = r + R~f, is equidistant 
relative to the surface on which the contact point is moving. 

In the case of potential forces, the moment Mp is expressed in terms of the potential U(r') = 
= U(r + Rry). This potential depends on the position of the ball's center of mass as follows Mp = 

or 

Remark 1 . In his treatise ^01 Routh obtained the equations of motion of the ball in semifixed axes and 
presented the cases when these equations can be analytically solved. Actually, in the majority of the subsequent 
publications E] Routh's results were just restated without any essential expansion. It should be noted that 
Routh was especially interested in the stability of particular solutions (e. g. a ball rotating around the vertical 
axis at the top of a surface of revolution). Here we are not giving the original form of Routh's equations. 
Equations (|2.3|) arc, in many respects, similar to the equations describing an arbitrary body motion on a plane 
and a sphere [S]- This allows us to consider many problems (e.g. those of integrability) from the single point 
of view. 

The integrals of motion. In the case of potential field with potential U{r + R*y), equa- 
tions Q2.3|) possess the integral of energy and the geometric integral 

H=^(M,u) + U(r + R'r), F 1=7 2 = l. (2.5) 

In the case of an arbitrary surface F(r) = 0, apart from these integrals, the system 1)2. 3 Jl has neither 
measure, nor two additional integrals, which are necessary for the integrability according to the last 
multiplier theory (the Euler-Jacobi theory). The system's behavior is chaotic. As it is shown later, 
in some cases there exists a measure and only one additional integral. In such a case chaos becomes 
"weaker". As it was noted by Routh, for a surface of revolution there exist two additional integrals, 
the system is integrable, and its behavior is regular. The reduced system becomes a Hamiltonian one 
after an appropriate change of time. 

Rolling on a surface of the second order. We are going now to revisit equations (|2,3|) for 
the case when the ball's center of mass is moving on a surface of the second order. The surface is 
defined by the equation 

(r + R<y,B-\r + R<y)) = 1, B = diag(&i, b 2 , 63) (2.6) 

(for an ellipsoid, b% > and are equal to the squares of the main semi-axes). Solving ()2.6|) for the 
position vector r gives 

Bt 

r + ff 7 = — 7 . (2.7) 
We obtain the equations of motion in terms of the variables M,y. 



^ = ~7T+S (M '^ 7 ' 7= ^T^ 7X bxB-^TxM)) (2.8) 

A ball on a rotating surface. Let us also consider the motion a ball on a surface rotating 
with constant angular velocity fi. One particular case of this problem (a plane and a sphere) was 
also investigated by Routh By analogy with the previous case, replacing the nonholonomic 

relation lj!T2J) by 

v + uxa = Qxr, (2.9) 

we get 

M = m(a x (u) x a) + a x (ft x r)) + Mp 

y v ' \ i) f £ 2 _ 10 ^ 

r + i?7 = u x Rj + fl x r. 

Here a = —R*y. 
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It should be noted that in the case of an arbitrary surface the system of six equations (|2. 10|) is 
not closed because the position vector of the point on the surface, r, is not expressed in terms of 7 
only; one should also introduce the equation for the angle of rotation of the surface around the fixed 
axis. Nevertheless, if the surface is axisymmetric, and the symmetry axis coincides with the axis of 
rotation, equations (|2.10|) get closed. We assume this to be fulfilled in further text. 

Equations (|2.2j) . 1)2.10(1 are in many respects similar to the equations, defining the rigid body 
rolling on a plane or a sphere, which we thoroughly investigated in the paper We shall use this 
fact while transferring the corresponding results onto systems (|2.2[) . (|2,l(Jj) . 



3. Ball's motion on a surface of revolution 

First of all, let us consider the cases of integrability of equations (|2.3fl . (|2.1()j) . associated with the 
rotational symmetry of the surface on which the ball rolls. Here we are using the technique introduced 
in [H] and concerned with the analysis of a certain reduced system in terms of new variables K\ , K 2 , 
-^3; 73- We shall also assume the surface to be rotating around the axis of symmetry with constant 
angular velocity fi = (0,0,0), £1 7^ 0. The particular case, $7 = 0, was discussed by Routh, who 
obtained the majority of the following results (although he missed some cases of equal interest). 

The equation of a surface of revolution with respect to an immovable reference frame may be 
written as 

H = (/(73) - Rhi, r 2 = (/(7 3 ) - #)72, 

U !- 73 2 , \ M 

r 3 = / (/(73) / (73)Jd73 - #73, 

where 7(73) is a certain function specifying the surface parametrization. The parametrization (|3.1|) is 
to be so chosen as to make the form of the reduced system as simple as possible. 

In the case being considered equations (|2.3j) . (|2. 10|) allow an invariant measure with density 

P = (/(73)) 3 3(73), where 5 ( 73 ) = /( 73 ) - i-^i/'( 73 ). (3.2) 
Apart from the invariant measure, the equations also have a simple symmetry field 

i=i " h A- ih A +ii A-<^' <3 - 3) 

which is caused by the rotational symmetry. 

In our book [3] we make a frequent use of the fact that to obtain the simplest form of the reduced 
system in the presence of symmetries, one needs to choose the most relevant integrals of the field of 
symmetries which, in this case, may be written as 

v our \t( \ /xM 3 + D(M,7)73 M ll2 - M 27l 

K x = (M,7)/(7 3 ), K 2 = ntu 3 = — — , K 3 = — , 73. (3.4) 

H + Lf Vl - 73 

In terms of the chosen variables, the reduced system takes the form 

73 = kK 3 , 

^ = kK ^{) K ^< l -i)^ (3.5) 

Ks = -k + ^ " 73/(M + D)K 2 ) + 

^(1-73)/ V /(I -7a) 

+ (l - |) ((// + 2D 7 3 2 )^ - 73 /(/x + 2D)K 2 ))^ 
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Ry/T 



. It is easy to show that this system of equations possesses an invariant measure 



where k 

+ D)g(j 3 ) 

with density p = k~ l . The system 1)3.5)1 can be explicitly integrated in the following way. 

Let us divide the second and the third equation of system (|3.5j) by 73. Then we obtain the 
nonautonomous system of two linear equations with the independent variable 73 



dKi 

c?7 3 



L 

73 



K 2 + 1 



I 
R 



dK 2 

c?7 3 



D 



73(1 



R 



fin). 



(3.6) 



This system of linear equations always possesses two integrals which are linear in K±, K 2 . The 
coefficients in the integrals are functions of 73, and in the general case cannot be obtained in the 
explicit (algebraic) form. Having divided the last equation of the system (j3.5|) by 73 and substituting 
into this equation the known solution of the system we obtain the explicit quadrature for ^3(73). 
Using the first equation from (|3.5|) we can obtain the expression for 73 (t). 
In the case O = 0, the system (|3.5|) has the integral of energy 



H = (M,w) 



1 



A/ 2 



+ 



(ji + D^Kx-fKz) 



^/ 2 (1 



ll) 



fi + D 



(3.7) 



The quadrature for 73 (t) may be obtained upon substitution of K 3 = k 1a / 3 into 1)3. 7 Jl . 
For the reduced system ()3.5)) the following theorem is valid: 

Theorem 1. For £1 = by the change of time kdt = dr, the system 1)3.5)1 can be represented in 
the Hamiltonian form 

(1 T ' 

^ = { Xl ,H}, x = ( l3 ,K l ,K 2 ,K 3 ). (3.8) 
The bracket is degenerate and specified by the relations 

{ l3 ,K 3 } = p + D, {K 1 ,K 3 } = (jjL + D)^K 2i {K 2 ,K 3 } = j r K 1 



(3.9) 



(for the other pairs of variables the brackets are zero). 

Proof is a straightforward exercise consisting in derivation of ()3.5)) (O = 0) from 1)3.8)1 and 
verification of the Jacobi identity for the bracket ()3.9)) . ■ 

One can assert that the equations for the linear integrals of the system H3.fi)) exactly coincide with 
the equations for the Casimir function of the bracket ()3.9)) . 

It should also be noted that for f2 = the system 1)3.5)) has "skew-symmetrical notation" , similar 
to the Hamiltonian form 1)3.8)1 . (|3.9)) 



d x i _ t ( \ OH , 



dr 



J 



(3.10) 



The matrix is of a somewhat more general form than that of the corresponding bracket ()3.9)) 



/ 









V 



u 



jp + D) 2 (73^1-/^2 
^ 2 /(l- 7 3 2 )^3 






-A 

((x + D)tlC 2 



Xu 




A 


jK x + Xv 



p + D 
J' 



\ 



(p + D)^K 2 + Xu 
^K x + Xv 



J 

(/x + D)((» + D^K, - /73GU + D)K 2 ) 

» 2 f 2 (i-ii)K 3 



(3.11) 
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Here A is an arbitrary function of 73, Ki. For A = 0, one again gets bracket (|3.9|) . and the tensor J\ 
satisfies the Jacobi identity. However, in the general case, J\ does not satisfy the Jacobi identity (it 
is not of the Poisson type). Nevertheless, if we assume 

A = ^(73)^3, (3.12) 

then the tensor A _1 Ja is a Poisson one, and the quantity A is a reducing multiplier (according to 
Chaplygin). 

Thus, we obtain the following Hamiltonian vector field 

v x = (A~ 1 J A )VF, 

for which divi> A 7^ (divAi^ = 0). The examples of Hamiltonian fields with nonzero divergency were 
almost not discussed earlier. One particular case of the Poisson tensor A -1 J\ was found by J. Hermans 
in |12| . Hermans used his own system of reduced variables which slightly differs from ours. 

To clarify the behavior of a ball on a surface of revolution, we will discuss the cases, when this 
surface is a paraboloid, a sphere, a cone, and a cylinder. These problems were investigated by Routh 
in ^01 for $7 = 0. Here we will expand the results for f! / 0. 

A paraboloid of revolution. Suppose the ball's center of mass is moving on the paraboloid 
of revolution z = c(x 2 + y 2 ). In equations (j3.1|) we, therefore, put 

7(73) = -^. (3.13) 

The density of the invariant measure 1)3. 2 j) (up to an unessential factor) can be written as 

p=\. (3.14) 

7s 

In this case, the two-dimensional system (|3.6j) reads 



£i = -!,(*- ±±£%n). ^ = -^K* 1 + <i±|^,n|. <s.i») 

^73 2c7fV 2oR7 3 / d 73 n + D\ 2cR^ > 

For the variable K 2 we obtain a homogeneous second-order linear equation whose coefficients are 
homogeneous functions of 73 (at $7^0) 



3, 



j,,, 1 r,, , D 1 u D ( 2 + ^73 ) n 
2 73 2-r /i + £) ^2 (^ + D)cR~ff 

Its general solution may be represented as 73 , a = const 
For the variable K\ from 1)3.15)1 we, similarly, obtain 

* - ^H 1 - + + ss^j^)- (3 - 17) 

Solving for the constants ci, C2 from (|3.16|) . (|3.17|) . we obtain the integrals of the system (|3.6j) . These 
integrals are linear in to K\, K 2 and have form 

^2 = , = 73 2c 73^1 + -fj ^2 + 73 



+ V ^73 V D73 2oR(2-i/)7|' 

77 -*Yo u , ^ + g r^ , of (/i + £>)(l-^) 2-1/ ' 

F 3 = ; — ^ -73 ( 2C73K1 + n ^ -RT 2 + ( 73 ^ , — 2 



+ 6 V ^73 " V'° Z) 73 2ci?(2 + 7/) 7 . 
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The product F 2 F% gives an algebraic quadratic integral. 

For $7 = the equations also have the integral of energy 

* = ^4^M* + *)' + ^ (3, 9 ) 

Remark 2. Some forms of surfaces (on which a ball is rolling) are investigated in the paper jS] (details 
are given below). In this paper it is shown that if a ball is rolling on a paraboloid of revolution, then the 
system (|2.3|) is reduced to a particular class of Fuchsian equations. Routh himself considered the case when not 
the contact point, but the center of mass is moving on a paraboloid. We should note again that in this case 
equations 1|2.3|) have algebraic integrals. The motion of a homogeneous ball on a surface of revolution was also 
studied by F. Noether [TT] . 

An axisymmetric ellipsoid. Consider the motion of a dynamically symmetrical ball when its 
center of mass is moving on a fixed ellipsoid (f2 = 0). In this case 

/(73) = Ttt, \ l r ( 3 - 2 °) 

- 7 3 2 ) + ^37 3 2 

where 61,62 are the squares of the ellipsoid principal semi-axes. The density of the invariant mea- 
sure ()3.2j) (up to a constant factor) is 

p= (6 1 (l- 7 |)+ 6 3 73 2 )" 3 . (3-21) 
In this case the variables Aq, JV2, K% are more convenient than the variables ()3.4|) . Here 

Aq = (M, 7 ), Ar 2 = _^_/( 7 3)( 7 3(M,7)-M 3 ), (3.22) 

which satisfy the system, similar to (|3.5|) for O = 

Ni = -kK 3 ^N 2 , N 2 = kKs-^-fNu 
73/ + D 

^3 = ~k 2 { ^ +D \ 9 , 2f3 N 2 (vf( l ~ 7s 2 )^i + (M + D) l3 N 2 ) , 
A» 73(1-73) / 



(3.23) 



where k 



{H + D)g 

Using (|3.20j) . we obtain two linear equations with independent variable 7 3 

^1= - {bl ~ h) N 2 ., <^ = ^ AT, (3.24) 

d ™ 6iV&i(l - 7 3 2 ) + &37 3 2 d ^ (v + DWh(l-ii) + b 3 ii 

It is easy to show that system (|3.24|) has a quadratic integral with constant coefficients 

F 2 = b i—Tn N ? + " b *) N l ( 3 - 25 ) 

Below we will show that this integral may be generalized to the case of a three-axial ellipsoid. 

The system (|3.24[) can be solved in terms of elementary functions. Depending on the sign of the 
difference 61 — 6 3 , the solution may be written as: 
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1. bi > 6 3 , a 2 = 



bi — b s 



> 1. 



Ni = a sin 99(73) + c 2 cos 1/2(73), N 2 = a \j ^ D (~ c i cos ^(73) + c 2 sin 1/3(73)) , 



^(73) = farctg 



13 



/x + £> 



(3.26) 



2. bx < b 3 , a 2 



bi 



63 — bi 



> 0. 



Nx = cit~ v + c 2 t u , N 2 = a< 



Ubi 



fi + D 



-ClT~ V + C 2 T V 



Tim) = 73 + \Ja 2 + 7|, 1/ 



(3.27) 



Here ci,c 2 = const. Solving for these constants, we can obtain linear integrals of motion. 

Historical comments. It is an interesting fact that neither Routh, nor his followers succeeded in ob- 
taining the simplest reduced equations (like l|3.24|l ) and solving the problem of a ball rolling on an ellipsoid 
of revolution in terms of elementary functions. To integrate the equations one must appropriately choose the 
reduced variables such as 




Fig. 2 



map 7 



A circular cone. In this case, due to the fact that the Gaussian 
VF 



is degenerate, one should use the components of the 



vector r (the position vector of the ball center) as the positional 
variables in equations (|2. 1U|) . For the cone 9 (Fig. |2J) we have 



73 = cos 9 = const , 71 



k 2 ri 



r 3 = kJ r'f + r%, k = t 



72 



k 2 r 2 



(3.28) 



The measure of equations (|2.3|) . (12.10(1 in which 7 is expressed in 
terms of r in accordance with (|3.28j) . can be written in explicit form: 



P 



k 2 



l-k 2 



Ro + r 3 



(3.29) 



For the reduced system let us choose the variables 

Dtt r 3 



07 = u 3 + 



c 2 = r 3 + 



k- 



VTTk 

In terms of these variables we obtain the equations 



R M 



2, ,2 



n ,7 



^l-n d<7 2 _ fT—j^ »Sl (r 3 2 



(3.30) 
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Prom these equations the following integrals can be easily obtained: 



\i + D \2R 



(3.31) 



A circular cylinder. The motion of a ball within a cylinder is a well-known problem which is 
usually used to illustrate some unrealistic conclusions, derived by means of nonholonomic mechanics. 
It can be shown that a homogeneous ball, moving within a vertical cylinder, on the average, is not 
rolling downwards due to the gravity force. Nevertheless, this physical fact may be observed while 
playing basketball, when the ball has almost hit the basket, but then rapidly jumps out of it, suddenly 
lifting upwards. The addition of a viscous friction to this nonholonomic system, which leads to a 
vertical drift, is analyzed in |Sj, where the explicit solution of this problem has also been obtained. 

~, ) , where R c is the cylinder radius (Fig. EJ), in 

Rc Rc ) 

terms of the variables (M, r) or (u>, r), the invariant measure's density is constant. 
The kinetic energy looks like 



For a cylinder, 7 



H = \{M, u) 



i(/i(cj, 7 ) 2 + + D)(f4 + (^172 - ^27i) 2 ))- 




The reduced system in terms of the variables a\ = 0J3, 02 = (w, 7) looks like 

Ruj 3 - R c n 



a 1 =oj 3 = 0, 



/ 

(To 



(Rc — R)R 



(3.32) 



Fig. 3 



Thus, we have two integrals 

= const, 



, \ Ruiz - RSI r 3 

(cj ' 7) — i — jh~R = const - 



(3.33) 



On writing u = (uj\, 0J2, ^§ — ), the second integral takes the form 

R r — R 



(u, r) = const, 



and the kinetic energy becomes 



2H = /i I (2, 7) + S 3 ^J + (/i + mR 2 )u 2 3 + (/i + mR 2 )^. 
Hence, the variable r 3 (responsible for the vertical displacement of the ball) can be easily found to be 

(2,7) 



^3 



/ 2g-Qi + m tf)Sg HiTl(Si| 



l_fJ_ 
fi + vtlR' 



it -to) 



where to is a constant depending on the initial conditions. It is clear that the average displacement of 
the ball equals zero even in the presence of the gravity field (see formula (|5.6|0 . 
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r. 







Let us now consider two different variants of problem concerning the 
rolling of a ball on a sphere. These problems were integrated and analyzed 
a by Routh [TO]. 

A ball on a rotating sphere. Consider a sphere of radius R s ro- 
tating around a certain axis with a constant angular velocity fl. Let R be 
the ball radius, a = — i?7, r = R s ~f (Fig. EJ). The equations of motion in 
the potential field with potential V(j) may be represented as 

Dl " = i^ft^' 7 ^ X 7 + R~TR { ^ 1){RW + RsQ) X 7 + 7 X %> 



Fig. 4 



7 



+ i? s S7) x 7 



(3.34) 

where D = mR 2 , D\ = fi + D, and // is the moment of inertia of the ball. 
Whether the ball rolls outside or inside the sphere is determined by the sign 
of R. Using the angular momentum vector M = D±u — 1)7(0;, 7), the first 
equation of (j3.34|) may be written as 



M 



x 7 )(n, 7 ) + 7(0,w x 7)) - fl,(0,7)(7 x «)) +7 x 



(3.35) 



Using u> = u> + -j^fi in (13.341) . we get 



7 



i? 



(3.36) 



a; x 7. 



This coincides with the original system ()3.34j) for fi = 0. Thus, it will be enough to consider only the 
case when f2 = 0. In terms of the variables (M,-f) equations ()3.34j) for Q = are written as 



M = 7 x 



dV 



7 



R 



(3.37) 



R s + R 



u x 7. 



They possess the integrals 



H=±(M,u,) + V('j), F\ = (M,7) 



The change of time t 
equation 7 



F 2 = 7^ = 1. (3.38) 
r transforms the second equation of ()3.37j) into an ordinary Poisson 



R s + R 

7x0?, and the potential is multiplied by some nonessential factor. Thus we set a 
system (|3.34|) for which we know the cases when it is integrable. For example, if V = ^(7, C7), 

C = diag(ci, C2, C3), one gets famous Neumann problem of the motion of a point on a sphere in the 
quadratic potential (equations (|3.34j) incorporate even more general situation, when (M,7) = c ^ 0, 
which corresponds to the Clebsch case jl]). 

The equations (|3.36j) . (|3.37j) for the motion of a ball on a sphere also manifests an analogy (dis- 
cussed in OE3) of problems, concerning the rolling of a homogeneous ball on a sphere in the gravity 
field, and the Lagrange case in the Euler- Poisson equations (the motion of a heavy dynamically sym- 
metrical top). Indeed, for a potential of the type V = ^(73), the systems (|3,3fi[) . (|3.37|) . because of 
the axial symmetry, have "Lagrangian integrals" M3 = const or W3 = const. 
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The rolling of a ball on a rotating sphere. Now suppose that the sphere on which a ball 
rolls rotates freely about its center. The dynamical equations can be written as: 

mi) = N, fid) = a x N, /x s f2 = — r x N, 
v + u>xa = Q,xr. 

Here u, fi, fi, fi s are the angular velocities and the moments of inertia of the ball and the sphere, 
respectively. Using the relations r = R s ~f, a = —R~f, for the quantities cD = — , 7, we can 
write 

(1 + D)u = Dj x (7 x u>), 7 = a) x 7, 



mg ( 3 - 39 ) 
A* 1 Ws ' 

By the change of time cZi — > a dt, a = const, this system can be reduced to equations (|3.36j) and, for 
this reason, is integrable. 

4. Rolling of a ball on surfaces of the second order 

Consider the dynamics of a ball in greater detail for the case when its center of mass is moving on a 
surface of the second order 

(r + Ry,B~\r + R 7 )) = 1, B = diag(&i, b 2 , b 3 ). (4.1) 

In this case the equations of motion are identical in form to (|2.8[) . 

It can be shown that these equations possess an invariant measure and a quadratic integral of the 

form 

p= (7,B 7 )" 2 , 

_ (7 x M,B- 1 (7 x M)) (4.2) 

F2 ~ ( 7 ,B 7 ) ' 

Here B is an arbitrary (nondegenerate) matrix. 

Comment I. The invariant measure's density was found rather easily, after the authors had obtained 
equations l|2.8|l for rolling of a homogeneous ball on an ellipsoid. These equations describe the Jacobi non- 
holonomic problem. The problem's name stems from the following fact. When the ball's radius is tending to 
zero it seems that we get the holonomic classical problem concerning the geodesic lines on an ellipsoid (solved 
by Jacobi in terms of elliptic functions). Apparently, one should get convinced in the correctness of such a 
limiting transition, which, however, does not prevent us from using the given terminology. One can treat the 
integral (|4.2|l as a generalization of the Ioachimstal quadratic integral in the Jacobi problem. Originally, the 
authors found this integral by numerical experiments, using the Poincare three-dimensional map in terms of the 
Andoyae-Deprit variables (L, G, H, I, g, h). These variables for nonholonomic mechanics were introduced in |Sj 
(see also earlier paper 

Fig. illustrates three-dimensional sections of a phase flow in the phase space (l,L/G,H,g) at a level of 
energy E = const. The secant plane is g = ir/2. It is seen that the levels F2 = const make the three-dimensional 
chaos "foliated" into two-dimensional chaotic surfaces. The fact that the motion on the two-dimensional surfaces 
F2 = const ensures that no additional integrals (necessary for intcgrability of the problem) exist. 

Let us consider various particular (may be degenerate) surfaces of the second order for which an 
integral of the type (|4.2|) exists. 

An elliptic (hyperbolic) paraboloid. Let the ball's center of mass move on an elliptic 
paraboloid defined by the equation 

f- + f- = 2z. (4.3) 
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Fig. 5. Examples of three-dimensional mappings (on the left) and the corresponding level surfaces of the integral 
F2 (on the right). All the mappings are constructed at a level of energy E = 1 for B = diag(l, 4, 9). To the 
frames on the left correspond the values F2 = 1.7, 2, 4, 8, 10 (from top to bottom). 
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Although in this case all the results may be obtained from the previous ones by means of a passage 
to the limit, we will derive them "of scratch". 



The Gaussian map (|2.7|) looks like: 



r 1 + R 7l = -b 1 ^, r 2 + R 72 = -b 2 ^, r 3 + R>y 3 = ^ ^ - , (4.4) 



and the equations motion (|2.8j) take the form 



M = --^(M, 7)7, 7 = ^2^ 7X ( 7 xB*( 7 xM)) ; 



fi + D" 1 n + D 

where B* = diag(6j" 1 , \ 0) is a degenerate matrix. 
The invariant measure density depends on 73 only 



P=\, (4-5) 

73 

and the quadratic integral (j4.2|) can be written as 

( 7 xM,B'( 7 xM)) 
-^2 = 5 • (4.6) 

73 

Comment II. The integrals 1)4. 2[) . 1)4.6)) . being quadratic with respect to the velocities (M or u;), depend 
on the positional variables in a rather complex way. Maybe for this reason the classics (in particular, Routh and 
F. Noether who obtained only particular results) did not find these integrals. As it has been already noted, the 
integrals 1)4. 2J) . 1)4.6)1 were originally found through numerical experiments. Their analytic form was obtained 
by means of the following considerations. 

As we have shown above, the problem concerning the rolling of a ball on a paraboloid of revolution such 
that bi = 62 is integrable and has two additional linear integrals. This integrals (in the absence of rotation 
fl = 0) may be written as follows 



h = 73^ 1 (< J l - ^273 + C273V4 - k), 

-vl-fc-i, f, r, , D ( 47 ) 

h=7 3 (01 - 0273 - 0273V 1 - k), k = ^ + D 

where o\ = W3, 02 = ( w j 7)- The product of these integrals is also an integral which is quadratic with respect 
to Oi , and a rational function of o\ , <Ji and 73 



(0273 _ 



2 



J = cr 2 (l - k) = — - 2a 2 ^- + ka 2 . (4.8) 

7 3 7s 73 

Let us eliminate, using the expression for energy 1)3. 7)1 . the term ka\ from the integral l)4.HJI . and consider the 
integral 

^2 = J + 2S=4-2^2^f+^ 2 . (4.9) 

7 3 73 

Substituting the expressions for 02 into 1)4. fl)) and isolating perfect squares, we obtain 

(72W3 - 73^2) 2 + (73^1 - 7l^ 3 ) 2 , . n , 

^2 = 7. ■ l 4 - 1 ^)] 

73 

The integral F2, written in such a form, is easily generalized to the case of an arbitrary paraboloid (61 ^ 62) 



T-(72^3 ~ 73W2) 2 + T-(73Wl ~ 71 w 3) 2 



F 2 = - g-f = 11 L_^Z ^ B^diagfcrS^.O 4.11 

73 73 

and any surface of the second order 1)4.20 . The integrals 1)4.2)1 . 1)4.6)1 may be used for stability analysis of 
stationary motions of a ball near the points of intersection of the surface with the principal axes. The symmetrical 
case of this problem was considered by Routh ^U] . 
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Motion of a ball on an elliptic cone. Suppose that the ball's center of mass is moving on 
the surface of an elliptic cone, defined by the equation 

(r c , B-Vc) = 0, B = diag(6i, 62,-1), (4.12) 

where r c = r + R'j are the coordinates of the center of mass, and 61, 62 are positive quantities such 
that \fb\ and \fb<i determine the slope of the generatrices with respect to the coordinate axes. Given 
the coordinates of the center of mass, we can calculate the normal to the surface 7 at this point as 
follows: 

^(B-Cu e) - (4i3) 

In our case r (or r c ) is not uniquely defined by 7 (because 7 is constant on a generatrix). Therefore, as 
phase variables, we will use not (M, 7) (as we did earlier), but (M, r c ). Upon substitution of 1)4.13)) 
into the equations of motion (|2.3[) . we get the equations of motion in terms of the variables M, r c : 



R 



(M x B~ l r r 



Oi + D)VP=^ B-ir c ) (4 U) 

M = ^ r (B _1 (Af x B-V c ), B -1 r c x (B" 1 ^ x M))B~ 1 r c . 

Equations 1)4.14(1 possess an energy integral 

1 , 2 DIM, B~ l r c f s 
H=— — -(M 2 + V 



2(^ + Z?) v MB-Vc, B-Ve)' 

and an invariant measure 

p= V(B- 1 r c , B-!r c ). 
Let us now make one more change of variables and time 

y = B~ 1 r c , dr = = dt. 

( M + i))V(B- 1 r c , B-Wc) 

This results in 

y> = B- 1 {M x y), 

M ' = 7 £ ^(B- X (M x y), y x (y x M))y 

which possesses two "natural" integrals: the energy integral 

2^ + Z)) 1 A* (y, y) ; ' 

and the geometrical integral 

(v, By) = o. 

The latter defines the surface in terms of the new variables along which the ball's center of mass is mov- 
ing. Moreover, equations (|4.15|) possess invariant measure with constant density. The generalization 
of the nontrivial integral 1)4.2)1 to the case of equations 1)4.15)1 looks like 

F 2 = ((Mxy),B- 1 (Mxy)), (4.16) 

or, in terms of the original physical variables, we set 

F 2 = ((Af x B _1 r c ), B _1 (M x B -1 r c )). 

The question of existence of one more additional integral for equations 1)4.15)1 remains open. 
Apparently, in the general case, when 61/62, it does not exist. 
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5. Motion of a ball on a cylindrical surface 

Let us consider the rolling of a ball whose center of mass is moving on a cylindrical surface. It can be 
shown that, in the absence of external fields, this system may be integrated by quadratures; but when 
an external force is directed along the cylinder generatrix, the equations are reduced to a Hamiltonian 
system with one and a half degree of freedom . 

Let us choose a fixed frame of reference with one axis (Oz) directed along the 
cylinder generatrix (see Fig. |f)J). In this normal vector is expressed as 



7= (71,72,0), 7i+7| = l. 



(5.1) 



Denote the projections of the normal to the center of mass and the position vector 
of the center of mass of the ball onto the normal cross-section by r = (n + R"fi , r2 + 
+ -R72), 7 = (71,72)- For these projections we have evident geometrical relations 

(r,j) = (7,7) = 0. 
Hence, we conclude that 7 is parallel to r, 

7 = A(7)f . 













V 


e 



Fig. 6 



The factor A(7) is completely determined by the geometry of the cylinder cross-section and does not 
depend on angular velocity. 

Using equations H2.3|) . we get the equations of motion for the ball on the cylindrical surface with 
the assumption that the cylinder is subject to a force (with potential U(z)) directed along a cylinder's 

generatrix (z = ^|): 
K 



M 



7 



M3 D 



dU 



M 3 



(5.2) 



Kl) e z x 7, 



(M x 7, e z ). 



Here not only the energy is conserved but also the projection of the (angular velocity) moment 
on the cylinder axis: 



H = \{M,u) + U{z), 



(5.3) 



F2 = M3 = (fi + D)uJ3 = const. 
Moreover, the system ()5.3|l possesses an invariant measure with density 

p(7) = A" 1 ( 7 ). (5.4) 

It follows from (|5.2|) that the equations for the vector 7 get uncoupled. Let us use parametrization 

71 = cos (p, 72 = sin ip. 

For the angle ip(t) we obtain the equation 



<P 



M3 
H + D 



A(cos <p, sin (p) = Q 1 (<f), 



(5.5) 



where Q((p) is, in the general case, a 27r-periodic function of <p, determined by the shape of a cylinder 
cross-section. 
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In the remaining equations of the system (|5.2|) we put 

K\ = M171 + M271, K 2 = M x 72 - M271, 

and replace the time (as an independent variable) by the angle ip, and thereby obtain the nonau- 
tonomous system with 27r-periodic coefficients 



J\ (5-6) 

dz Qm r, 

K2- 



dip fi + D 
This system has an integral of energy 

In the case of the gravity field, U (z) = mgz and equation (|5.6|) are integrated by quadratures: 

Ki((p) = —vmg J sin v{t — (p)Q{r)dr + vA cos vip + vB sin ucp, 

fo 

<p 

K2( ( p) = —mg J cos v(t — (p)Q(r)dT + A sin vtp — B cos u<p, 
2 _ M 



(5i 



where A 5 are constants and u 

fx + D 

Let us show that the integrals in ()5.8|) are bounded functions. We expand the function Q(r) in a 
Fourier series 

Q(r) = J2Qne inT . (5.9) 

nez 

The integrals in the expressions for i^i(y) and -^(v 9 ) i R (|5.8|) may be considered as the real and 
imaginary parts of the integral 

! e ivr Q(r) dr= fj2 Qne i{n+u)r dr. (5.10) 

Using the well known theorems of the Fourier analysis and the fact that n + v ^ (as far as < v < 1) 
we put the integral under the sum sign and integrate the series term-by-term: 



/ Y Q n e i{n+U)T dr = Y Qn A n+ ^ T . (5.11) 
/ i(n + u) 



It is evident that the series obtained converges to a certain quasiperiodic function, hence, K\{ip) and 
K2((f) are bounded. This fact and also the conservation of energy in the reduced system (|5.7[) ensure 
that z{ip) is bounded. 

Thus, when the ball is rolling on an absolutely rough cylindrical surface of an arbitrary cross- 
section in the gravity field, the vertical secular drift is not observed. 

The time dependence of the angle (and, hence, of all the other functions) is described by (|5.5|) . 
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An elliptic (hyperbolic) cylinder. Let us consider in greater detail a particular case, i.e., 
the ball's center of mass is moving on an elliptic cylinder with cross-section is defined by the equation 



h b 2 



(5.12) 



We have 



ri + Rji = — , r 2 + Rj 2 = = , ==, r 3 = Rz, 



V&i7i + hii 



therefore, 



H~r) = Rh f h 7)3/ \ B = diag(6 1 ,6 2 ,0), 
&1&2 

(H + D)bib2 



We should mention an important distinction existing between an elliptic and a circular cylinder 
(see above): in the case of an elliptic cylinder the dependency of dynamical variables Kx,K 2 ,z is 
defined by two frequencies ui = 1, u 2 = v, instead of a single frequency, as it happens in case of 
a circular cylinder. Thus, the integrals in (|5.8|) contain quasiperiodic functions; the integrals have 
very complicated nature, their analytical properties are thoroughly discussed in Graphs z((p) for 
various initial values of K\ and K2 are shown in Fig. [Tj The main result is that whatever the ratio 
of the frequencies, the quantities K\ and K2, and, therefore, the displacement z execute bounded, 
quasiperiodic oscillations. It is the main result of this construction. 

The authors express are deeply thankful to V. V. Kozlov for useful remarks and discussions. 




Fig. 7. The (^-dependence of the vertical coordinate of the point of contact z for various initial values K\, K2, z. 
For this figure the other parameters are: E = 1, fi = 1, D = l(v = 2 -1 / 2 ), bi = 1, 62 = 2, R = 1. 
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Table 1. The rolling of a ball on a surface 



surface type 


cylindrical surface 


surface of the second order 


surface of revolution 


ellipsoid, hyperboloid, 
paraboloid 


cone of the second order 


arbitrary 
surface 


ellipsoid, 
hyperboloid 


paraboloid, cone, cylinder 


sphere 


measure 


measure exists 


p = ( 7 ,B 7 )- 2 




P=(/(73)) 3 (/(7 3 )-^^/'(7a)) 


additional 
integrals 


system is integrable 
by quadratures 


(7xM,B- 1 ( 7 xM)) 

— const 

( 7 ,B 7 ) 

(one integral) 


((M x B-V.J^-^M x B"V)) = 
— const 
(one integral) 


two linear 
integrals, denned 
by the system of 
linear equations 


there exist two linear integrals that can be expressed in 
terms of elementary functions 


Harniltonianity 


nothing is known about the Harniltonianity of 
these systems 


upon change of time (prescribed by the reducing multiplier) 
the reduced system becomes Hamiltonian 


authors 


A.V.Borisov, I.S.Mamaev, A.A.Kilin (2001) 


E.Routh (1884) 


A.V.Borisov, 
I.S.Mamaev, 
A.A.Kihn 
(2001) 


E.Routh (1884) 


generalizations 
and remarks 


integrable addition 
of the gravity field 
along a cylinder 
generatrix is possible 




A.V.Borisov, I.S.Mamaev, and A.A.Kilin have shown the 
integrability in terms of elementary functions of the case of 
the ball rolling on the ellipsoid of revolution; have found an 
invariant measure for an arbitrary surface of revolution; for 
a paraboloid, a cone, and a cylinder they have shown the 
integrability of the case, when the ball is rolling on a 
surface rotating around its axis of symmetry 


the problems ol the 
rolling of a ball on 
an unconstrained 
and rotating sphere 

are also solved. 
System also allows 
integrable additions 
of potentials. 



Remark. The cases when the tensor invariants exist are indicated by gray color in the table. The partial filling corresponds to the uncomplete set of integrals. 



